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Abstract 

O, 

CN I We study a transmission coefficient of graphene nanoribbons with a top gate which acts as an 

^ ' obhque barrier. Using a Green function method based on the Dirac-hke equation, scattering among 
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transverse modes due to the obhque barrier is taken into account numericahy. In contrast to the 2- 
dimensional graphene sheet, we find that the pattern of transmission in graphene ribbons depends 



' strongly on the electronic structure in the region of the barrier. Consequently, irregular structures 



in the transmission coefficient are predicted while perfect transmission is still calculated in the case 
of metallic graphene independently of angle and length of the oblique barrier. 
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I. INTRODUCTION 



Recently, there has been considerable interest in transport on a graphene sheet, a single 



OS 



atomic layer usually pulled out of bulk graphite. [iHc 

closely-packed honeycomb structures, electrons in it behave like massless Dirac fermions. 
One of its fascinating properties is the so-called Klein tunneling where perfect penetration oc- 



Si. Due to its unique two-dimensional 

m 



curs independently of potential barrier height, contrast to the conventional, non-relativistic 
tunneling. {g-S] This relativistic effect is basically originated from the gapless electronic dis- 
persion, which in turn leads to the connection between electron and hole states in graphene. 
A sufficiently strong potential, being repulsive for electrons, is attractive for holes and gives 
rise to hole states in the barrier to form channels through which electrons can penetrate the 
barrier. [3] 

The idea of the Klein tunneling realized on graphene sheets was suggested by Katsnel- 
son et a/.joj] and then several experimental attempts have been made to demonstrate the 



perfect transmission. |10l-|l2| Up to now, the perfect transmission is addressed in terms of 



T = (1) 



sudden phase shift of conductance as a function of magnetic field. [12| However, a more 
direct evidence for Dirac particles may be the incident angle-dependence of transmission 
coefficient, jol which is not realized experimentally yet. For incident angle 6 and barrier 
length D the transmission coefficient T of Dirac particles is given by, 

cos^ 6 
1 - cos^{kD) sin^ 9 

showing oscillating behavior as a function of incident momentum k, angle, and barrier length. 
However, in the case of a 2-dimensional graphene sheet it is hard to adjust the incident angle 
6 because there are randomly directed particles which result in averaging Eq. ([T]) over the 
incident angle. A collimation-gate method to remove the randomness of electrons may be 
not appropriate because it requires proximate implementation to a potential barrier less 
than a mean free path. 

A good candidate to resolve the problem is graphene nanoribbons (GNRs) because elec- 
trons in each transverse mode propagate along its axis (i.e., 6 = 0). Then, we can adjust the 
incident angle 9 definitely by adopting an oblique potential barrier with respect to the axis 
of GNRs. However, in this case the transverse momentum is quantized and correspondingly 
Dirac particles become massive depending on their occupation to transverse modes. So, it 
is interesting to ask about what is the transmission coefficient of GNRs with an oblique 



potential barrier and to see whether the perfect transmission occurs in ribbon structures 
even under the obhque potential barrier. 

In this work, the Klein tunneling is investigated in graphene nanoribbons (GNRs) with 
oblique barriers theoretically. Based on Green functions taking account scattering among 
transverse modes we show that the transmission coefficient also oscillates, however, with rich 
structures as a function of incident energy, angle, and barrier length. Through the analysis 
of local density of states, it is found that rich structures in the transmission coefficients are 
resulted from detailed hole states formed in the region of the potential barrier and associated 
interference. We also examine effects of inelastic scattering for the experimental realization 
of the Kleining tunneling by employing the level broadening in the simplest approximation. 



II. MODEL AND METHOD 



To examine the Klein tunneling we assume a mono-layered graphene ribbon along y- 
direction with armchair boundaries as shown in Fig. [TJ It is well known that a graphene 
ribbon can be metal or semiconductor depending on its width W. 13| We set W = Nao/2 
where 2N is a total number of carbon atoms in a unit cell and ao = 2.46A is the graphene 
lattice constant. Then, the system becomes a metal if + 1 is a multiple of three, otherwise 
semiconductors. 

The system under consideration is composed of two electrodes and a central device part. 
The device part has a finite graphene ribbon with a gate electrode while the lower and upper 
electrodes are assumed to be semi-infinite perfect GNRs. We describe the system with the 
4x4 Dirac-like equation. This effective-mass equation is known to give accurate low-energy 



properties of graphene. 
reads as, 



131 ] With a potential V{x,y) induced by a gate, the Hamiltonian 



^xPx ^yVy 





+ lV{x,y) 



Px ^yVy 



E 

ip = — ij; 

Vf 



(2) 



where a^^y is a Pauli matrix, vp ~ lO^m/sec the Fermi velocity, p^. = —ikd^ and Py = —ihdy, 
respectively. When a gate electrode is tilted with an angle 6 as shown in Fig. [H the 
potential barrier is generally a function of x— and j/— coordinates and can be modeled as 
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FIG. 1: A schematic drawing of the device is shown. The obhque gate with an angle 9 is assumed 
in a graphene nanoribbon having armchair boundaries. D and W represent length and width of 
the oblique potential barrier, respectively. 



V{x,y) = V{y - {x-W/2)ta,n6) where 

/ 

V{y) = Fo 



\^|(l-sin7r 



for \y\<\{D-d) 
for \y\>\{D + d) 
otherwise 



(3) 



Here, we introduce a transition region of the length d to reduce numerical error of a finite 
difference method used in the foUowings as well as for the effective mass equation of Eq. ([2]) 
to be valid. 

In the absence of a potential barrier V{x^y), since the system is a perfect GNR one can 



solve the equation of motion ana 
conditions have been formulated 



Really. For instance, the appropriate armchair boundary 
13| and used to examine bound states [l^. Wavefunctions 



are known to be plane waves along both x- and y directions as. 



00 

2 



(4) 



where Znk = VQu — ik/ \/qn + ik and 7 = ±1 denotes conduction(+) and valence (—) bands 
of graphene, respectively. The boundary conditions yield the following quantization for a 
wave vector in x-direction; 



2tt 



n 



n 



integer 



(5) 



and electronic energy is given by E = 7^fFV'^l~+~^ for a propagating wave in ?/-direction 
with its wave vector k. 

In general an incident wave with a certain transverse mode from the lower to upper leads 
is scattered to other modes due to a perturbing potential in the device region. Consequently, 
one should take into account multiple transverse modes to resolve a scattering problem. In 
our case, we choose the following basis function to describe scattered waves 



oiq-nX 



-76 



-iqnX 



\ 



-iq„x 



(6) 



and adopt Green function approach to calculate the transmission coefficient. Actually, Eq. 
dH]) denotes eigenf unctions of an infinite GNR, namely Eq. (jlj) at A; = 0. So, since they 
are orthonormal to each other, i.e., {(pl \ (jP^,) = 5nn'5-n' ^^"^ satisfy the armchair boundary 
conditions, we can expand the wavefunction of Eq. ([2]) as ip{x,y) = J2n'yX"-i(y)'^ni.^)- 
Resulting equation of motion for Xn-yiv) is summarized as, in a matrix form. 



El + -D{y)+h-^]x{y) = 



(7) 



where xiu) is a column matrix with its component Xn-~i{y) and matrices D and b are given 

by, 

7+y 



bn7,n'7' 



2 

y - 7 



(8) 



Next, we implement a lattice version of Eq. ([7]) by replacing d/dy with a finite difference 
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on a uniform grid. Detailed results are shown as, especially in a block tridiagonal form, 
/ \ / • \ 



-b a„_i b 
• -b a™ 
■ ■ -b a, 




b 

■m+l 



■ ■ 
b ■ 



V 



x{ym-i) 

XiVm) 
XiVm+l] 



J \ 







(9) 



where a^ = El + D{ym) and b = b/2A with a grid spacing A. Following a standard Green 
function technique, we define a retarded Green function of the system as inverse of the 
matrix in Eq. (Q with slightly shifted energy E E + it]. Especially, we are interested in 
the Green function of the device part which is represented by grid points m = 0,...,M — 1. 
For the Green function of the device we truncate the lower- and upper-lead parts from the 
original matrix. (l5] Then effects of the truncation manifest itself to self-energies as. 



/ ao-Si b ■ 
— b ai b 



G = 1. 



(10) 



• -b slm-2 b 
V ■ ■ ■ -b a^-i-Sc; / 

Here, S^'^ is called surface self-energy from the lower and upper leads, respectively, and 
is related to its Green function, Gl,u through T^l^u = —bG^, [/(rric, mc)b where rric is the 
adjacent index to the device part. Due to the block tridiagonal form, the Green function for 
each leads satisfies the quadratic matrix equation, [a + hG L,u{^c,^c)^]GL,u{^c,mc) = 1. 
Then, using Eq. ([8]) we obtain. 



A\Eyh^vl-ql 



{E + -fhVFqn)6nn'S-yY. 



(11) 



From the calculated Green function, the local density of states (LDOS) at the index m 
can be found by 



DOS{m, E) = ImTrG(m, m 

TT 



(12) 



where G(m, m) is a block matrix of the diagonal Green function at the index m. And the low- 
field conductance a of the graphene ribbon can be calculated using the Landauer-Biittiker 



formula 16 1 

2e2 /"^ / df 



h /_V dE 



T{E)dE (13) 



where / = {1 + exp{{E — fi)/kBT})~^ is the Fermi-Dirac distribution function with the 
chemical potential /i and T{E) is the transmission coefficient. In terms of the most upper- 
right component of Green function G(0,M — 1) from Eq. ffTOj) the transmission coefficient 
can be expressed as, 



T{E) = Tr[G(0,M-l)r^G^(0,M-l)r^] (14) 



where T^-^ = i[^L,u - 



L,Ui- 



III. NUMERICAL RESULTS 



In this section, we numerically illustrate solutions of the Green function representing 
scattered waves from an oblique potential barrier and related transport properties. We 
consider a typical size of the device structure which may be realized experimentally; for 
instance, W = 99ao, D = 60ao, and d = SOao- To include a pattern of an oblique potential 
barrier, large length of the device part is chosen. Actually since the system is in equilibrium, 
the transmission coefficients are independent of the total simulation length as long as the 
scattering potential is described properly in it. We use 100 transverse modes for the accurate 
description of scattering waves and a grid spacing of A = 2A. So, for the total simulation 
length of 260ao a dimension of the matrix Eq. in ( ITOl) is about 32000. Instead of a full- 
storage scheme we solve the matrix equation with a standard tridiagonal inversion which is 
very efficient in computational demands. 

In Fig. m we show calculated transmission coefficients through potential barriers at 
different oblique angles and compare results for the cases of metallic and semiconducting 
graphene ribbons. For a perpendicular barrier {6 = 0) to incident waves, one can find the 
smooth variation of the transmission coefficient as a function of incident energy and sudden 
rising at every new occupation to a transverse mode. In the case of the perfect transmission, 
the staircase patterns are expected and thus the deviated ones in the figures imply back- 
scattered electrons due to the potential barrier. It is found that the perfect transmission 
occurs only in the lowest subband in the metallic case. Actually, this subband has Qn = 
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FIG. 2: Calculated transmission coefficients as a function of incident energy are compared for (a) 
semiconducting {W = 99o) and (b) metallic iW = 98ao) cases at oblique angles 9 = 0, 7r/12, and 
7r/4, respectively. We choose D = 60ao, d = 30ao, Vq = 0.5eV for potential barrier parameters 
with a grid spacing A = 2A. 



and thereby has the hnear dispersion hke that in a 2-dimensional case, which results in the 
perfect penetration as imphed in Eq. ([T]). 

As the obhque angle 6 increases, rich structures are found as shown in the second and 
third rows of Fig. [2l Comparing results at ^ = 0, one can see that the transmission coefficient 
is very sensitive to the incident energy; many peaks and dips appear in the small range of 
energy. For ^ = 7r/4 cases, calculation results are even similar to resonant tunneling in usual 
tunneling problems. On the other hand, it is very interesting to observe robust behavior 
of the lowest subband for the metallic case against the obliqueness of the potential barrier. 
The perfect transmission is retained within the presence of a single mode, i.e., below the 
onset of the second transverse mode regardless of the oblique angle. 

This behavior can be understood when we consider the scattering problem in terms of the 
Fermi golden rule, the first order perturbation based on eigenstates of Eq. for infinite 
GNRs. If the incident energy is small to occupy only the lowest mode, electrons cannot be 
scattered by the oblique barrier to higher modes due to the energy conservation. This in 
turn means that for an incident momentum k electrons can be scattered forwardly {k) or 
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FIG. 3: For semiconducting GNR iyV = 99o) with the obhque potential barrier {0 = '7r/4), we plot 
calculated transmission coefficients in (a) and local density of states in (6), respectively, for various 
barrier lengths. All curves are vertically offset for clarity and other parameters are the same as 
those in Fig. [2j 

backwardly {—k) within the same mode. However, in the case of a metallic graphene ribbon 
the lowest mode has a momentum of g„ = and thus z^n in Eq- (H]) has exactly the out-of 
phase for the backward-scattering relative to the incident one. Consequently, there is no the 
backward scattering in the metallic GNRs independently of the oblique angle. This is not 
the case for semiconducting GNRs because the lowest mode has gn 7^ and the out-of phase 
cannot be achieved for backward-scattered waves. 

The barrier-length dependence of the transmission coefficients is examined in Fig. |3]-(a). 
Contrast to oscillating behavior of the 2-dimensional graphene sheet as in Eq. ([1]), more 
peaks and dips appear in the transmission coefficient and are found to be blue-shifted as we 
increase the barrier length. For understanding of the calculated results we first note the blue 
shifted behavior. Interference effects are not appropriate to explain it because wavelengths 
of incident waves usually obey the geometrical relation of kD = — q^D = constant, 

which implies the red-shifted patterns with increase of the barrier length. So, as an usual 
explanation of the Klein tunneling we think about the connection to hole states in the 
barrier region. For this we examine the local density of states at the midst of the barrier 
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FIG. 4: For semiconducting GNR (W = 99ao) with the obhque potential barrier {9 = vr/4), gray- 
scaled local density of states are depicted as a function of energy along the GNR axis. D = 60ao, 
d = 30ao, and Vq = 0.5eV are chosen. 

and show calculated results in Fig. [3]-(b). Interestingly one can find that the variation 
of the local density of states is similar to that of the transmission as a whole. That is, 
peaks representing localized states become blue-shifted and more peaks at lower energy 
side appear as the barrier length increases. This is in accordance with usual behavior of 
states in a quantum well as its size is varied. Furthermore one can see that the interval 
between peaks is approximately equal to the energy- level difference in a quantum well, i.e., 
AE = hvF/ma:x{D,W} inferred from Eq. ([5]). Consequently, we attribute the variation 
of the transmission coefficient to that of hole state in the barrier region and peaks in the 
density of states contribute largely to the transmission coefficient. 

However, it is noted that all peaks in the density of states are not reflected in the trans- 
mission. In Fig. m we show the local density of states as a function of energy along the 
device axis. One can see that high density of states in the barrier region causes a peak in the 
transmission coefficient. On the other hand, the dense region in the local density of states, 
for instance around E = 0.03eV, corresponds to a dip in the transmission, which means 
that detailed transmission is also affected by interference; We attribute this to destructive 
interference. 

Now, we examine effects of inelastic scattering on the transmission occurring in graphene 
ribbons possibly from phonon, edge roughness, and impurities. To estimate its effects 
roughly, we adopt the simplest approximation where a diagonal self-energy representing 
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FIG. 5: Effects of the level broadening are examined for metallic GNRs with an oblique barrier 9 = 
7r/4 for 7/ = (solid line), O.lmeV (dotted), ImeV (dashed), and lOmeV (dot-dahsed), respectively. 
We use different sizes of the barrier in each panel while other calculation parameters are the same 
as those in Fig. [2j 



scattering is added to the Hamiltonian independently of position and energy. This effect is 
equivalent to the level broadening by setting a finite value ii in the retarded Green function. 

nn 

The value t] is related to the scattering rate 7 via the relation of 1] = ^7.p/7|, 118| According 
to Ref. 19, the scattering rates 7 are calculated to have values ranging about 1 ~ lOOTHz, 
equal to rj = 0.66 ~ 66meV at T = 300K for W = 5nm or smaller values for wider ribbons 
and lower temperature. 

In Fig. [5] we show effects of the level broadening on the transmission coefficients for 
metallic graphene ribbons at = 7r/4. It is found that the pattern of transmission becomes 
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smeared with increasing the broadening together with accompanying many new dips in the 
region of the perfect transmission, and eventually has no structures up to r/ = lOmeV. We 
also find similar behavior for different size of GNRs as shown in the second and third rows of 
Fig. [21 implying that the new dips are irrelevant of the energy difference between transverse 
modes. On the other hand, in the case of = 0, no new peaks are developed while the 
transmission coefficient are still suppressed with increase of the broadening (not shown in 
the figure). Through the analysis of the local density of states, we see that positions of the 
dips correspond to those of abundant density of states in the barrier region. 

According to the Fermi golden rule, the level broadening gives rise to scattering to 
more diverse states by releasing the energy conservation from the delta function to the 
Lorentzian one. So, the scattering matrix is averaged over the energy interval of ~ 
min{?7, hvp/D, hvp/WtanO}. For this reason, the out-of phase of the back-scattered waves 
is no longer achieved and thereby the transmission coefficient becomes suppressed. How- 
ever, we find that the newly developed dips reflecting more strong suppression cannot be 
explained by the Fermi golden rule and may be attributed to higher-order perturbations 
representing multiple scattering and associated interference. As a results, since the pattern 
of transmission is washed out for the level broadening about rj = lOmeV, it is necessary to 
reduce scattering rates for the experimental realization of Dirac particles, for instance, by 
lowering temperature and preparing clean GNRs. 



IV. SUMMARY 



In this work, we studied the transmission coefficient of graphene ribbons with an oblique 
barrier based on the Dirac-like equation. In contrast to the 2-dimensional graphene sheet, 
the transmission in graphene ribbons is found to depend strongly on the electronic structure 
in the region of barriers. Consequently, irregular structures in the transmission coefficient 
are predicted, however the perfect transmission is shown in the case of metallic graphene 
independently of angle, width, length of oblique barriers. This behavior still demonstrates 
the nature of Dirac particles in graphene ribbons because very large potential barrier is 
assumed. We also examine effects of scattering by employing the simplest approximation. 
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